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^>stract 



We^aoint out that nontrivial quantumstatistics of vortices in pla- 
nEjE^uperfluid systems has its origin in the Schrodinger equation 
raCS^r then the Chern-Simons term. Vortices in superfluid he- 
films are not anyons because helium films are compressible. 
C^^ntum Hall effect layers admit vortices with nontrivial statis- 
tical interactions thanks to their incompressibility. We consider 
type 11 planar superconductors described by time-dependent 
noslmear Schrodinger equation. Inside of a vortex core super- 
conductor is in the normal state. Because of that vortices in 
superconducting films are anyons. We consider also vortex loops 
itf^dimensions but a phase factor due to linking of vortex loops' 
worldshccts is not a topological invariant. 
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"-^Jteractions of vortices in superfluids described by the non- 
linear Schrodinger equation are by now a well studied problem 
jl], ||, ||, |[ ^| . What has not been studied in detail is the quantum 
th^tjry of vortices. In jij there has been proposed a topological 
Latfrdau-Ginzburg theory for vortices in superfluid helium. The 
theory is equivalent to the standard nonlinear Schrodinger equa- 
tugnj provided the topology of vortex lines is determined. The 
authors of Q have derived, in 3- dimensional case. Berry phase 
picked up by the microscopic wave-function due to vortex-loops' 
motion. 

In this paper we consider mainly planar systems which may 
be relevant to description of thin superfluid films. It is shown 
that there is a Berry phase proportial to the area enclosed by 
the vortex trajectory. Mutual interactions are considered and 
it appears that vortices might be anyons if a superfluid helium 
film were not compressible. Thus we try to find a model with 
an incompressible condensate. The phenomenological models for 
quantum Hall effect appear to be a well known example. How- 
ever we argue that Chern-Simons term is not essential for trans- 
mutation of vortices' quantumstatistics. In fact the statistical 
interactions come solely from the Schrodinger term. Finally we 
consider type II superconductors. Thanks to the normal current 
carriers inside a vortex core a defect of superfluid charge inside 
a vortex with respect to uniform background can be nonzero. It 
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is shown that vortices in superconducting films are anyons. 

Vortices in superfluid helium films. In the local limit 
the Lagrangian density of the model for interacting hard-core 
bosons is B 
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where VP is the condensate wave-function defined on the plane. 
We assume quantisation of the model by path integrals does 
make sense. It is convenient to rescale the fields and coordinates 
in the above formula 
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and subsequently multiply the whole Lagrangian density by 
to obtain the dimensionlcss Lagrangian 



A field equation is 
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id t ij) = -v 2 ip + - i)ip 
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The model admits vortex solutions of the form ip — 
h n (r) exp inQ. n is the winding number. The profile function sat- 
ifies the equation 
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Asymptotices close to the origin and at infinity are respectively 
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Energy of the vortex is regular at the origin but logarithmi- 
cally divergent at infinity. The divergence can be regularised by 
putting the system into a finite-size box. Energy of a vortex- 
ant ivort ex pair is also finite since the topology of such a config- 
uration is trivial - there is no slowly falling down gradient of the 
phase of the scalar field at infinity. 

Effective Lagrangian for vortices. Let me first describe 
derivation of an effective Lagrangian for widely separated vor- 
tices. The field configuration can be well approximated by the 
product Ansatz 



ip = Y[h[x - X v (t)]exp[ia v Q(x - X v (t))] 



(7) 
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where the product runs over particular vortices. a v is a sign 
of circulation of the v-th vortex and 9(x) is the azimutal angle 
around zero of its argument. 

The Lagrangian density (Ef) can be rewritten as 



L 
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where we have replaced -0 = P 2 expi%. With the use of the 
product Ansatz (0) the Lagrangian takes the form 



L = -Eq- J d 2 x[ P dtx + p 
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The first term is minus the net energy due to nonzero winding 
number and to core contributions. The second term contains 
information about vortex interactions. This term is regularised 
by p - the modulus of the scalar field. The modulus is close to 1 
almost everywhere except the cores of vortices. Loosely speaking 
the phase gradients and the time derivative of the phase in the 
second term have to be integrated over the whole plane except 
vortex cores. The first term to be integrated can be rewritten as 



J d 2 x pdtQ[x — x v (t)] 
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If there were only one vortex a careful treatement of the inte- 
gration of the multivalued phase would lead to S 



L\" 



x v X v (t) x X v (t) 



(11) 



Upon quantisation by path integrals this factor would lead to 
a phase picked up by the wave-function. If a vortex trajectory 
were closed than 
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2na v x " area enclosed by trajectory" 
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The area should be taken with a positive sign for a clockwise 
motion and with a negative sign for an anticlockwise motion. 
For trajectories like, say, the numeral "8" one part of the con- 
tour gives a positive contribution and the other one negative. An 
equivalent interpretation is that the phase picked up by the wave 
function is proportional to the total mass of the superfiuid en- 
closed by the contour. This is the whole story in the one-vortex 
case. 

If there were more than one vortex an additional term would 
arise from (RG) in addition to (O) 
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where 



M = 2tt / rdr [1 - hj(r 
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is the defect of the superfiuid mass with respect to the uniform 
background. With the asymptotics (^) it becomes clear that the 
integral is logarythmically divergent. Thus there is no dilute vor- 
tex gas limit in superfiuid helium in which vortices might be 
anyons. 



Discussion. There is no long-range statistical interaction be- 
tween vortices in the nonlinear Schrodinger equation. Compress- 
ibility of superfiuid manifests itself in divergence of the integral 
in Eq. (|lj) . In the incompressible case this integral would be con- 
vergent but, as it was shown in [Q, two-dimensional superfiuid 
films are compressible. It is so because of the contribution to 
the energy density from the term | Vx | 2 which is nontrivial in 
the case of nonzero net superfluidity. In other words, the term 
^rh n in Eq.(|J) gives rise to the long-range term in the second of 
Eqs @. This tcr m might be matched if the local interaction in 
the model ([!]) were replaced by a two-body repulsive potential 
V(R) = -S^ with some critical value of the constant g. Such a 
model would have nothing to do with superfiuid helium as long- 
range interactions are attractive there. However it shows that 
one can have anyons in a model just because of the Schrodinger 
term and without any explicit parity-breaking terms like Chcrn- 
Simons interaction. The Schrodinger term has a remarkable form 
very similar to the Berry phase in Quantum Mechanics. If we 
consider evolution of the condensate wave-function due to slow 
motion of a vortex the contribution of the Schrodinger term to 
the effective action is 



r * 2 '-> d 
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where tp is the condensate wave-function dependent on a vortex 
position X. <> is the standard scalar product. 

Thus an obstacle is identified. The core must be well-defined 
and for the core to be well-defined we have to make the fluid less 
compressible. The contribution from the phase-gradient energy 
can be matched by gauging the model d^x ~ * ®nX + iA^. An 
example are the models with Chern-Simons interaction being 
a phenomenological description of the quantum Hall effect. In 
these nonrelativistic models jll], |l3|, [li], |l5| the scalar field of 
the nonlinear Schrodinger equation is minimally coupled to the 
C-S field 

d^ip -> - ia^ip . (16) 

The kinetics of the gauge field is governed by the Chern-Simons 
term 
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All the cited models possess the Bogomolny limit |l(| and in this 
limit there are static multivortex solutions. 

The dynamics of CS vortices in the Bogomolny limit has been 
intensively studied in relativistic models ^ |l^| and also in the 
nonrelativistic context Jl2|, |l5|j . One starts from a static multi- 
vortex solution in the Coulomb gauge. In this gauge the phase 
of the scalar field is a sum of azimutal angles exactly the same 
as in (0). Since the gauge field is transverse the CS term does 
not contribute to the effective Lagrangian. The only contribu- 
tion comes from the Schrodinger term iifi*dtip- This is exactly 
the same as in the nonlinear Schrodinger equation. A different 
choice of the gauge would not change the topology of the phase 
of the scalar field - the nontrivial phase factor for a closed tra- 
jectory would still come solely from the Schrodinger term. Thus 
the nontrivial statistics of vortices comes from the Schrodinger 
term and not from the CS term. The CS term transmutes the 
quantumstatistics of the elementary field quanta only. The sta- 
tistical interaction between vortices arises because there is a def- 
inite mass defect inside the closed trajectory of a vortex (or a 
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mass excess for Jackiw-Pi solitons). The Gauss' law establishes 
connection between matter density and magnetic field 



kB = i/)*ip 
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The mass defect is quantised because the magnetic flux of the 
vortex is quantised. The interaction term at long distances is up 
to a numerical factor proportional to 
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The prefactor ft is the same as in the CS term (|17]). This coin- 
cidence might suggest that the CS term is responsible for the 
statistical interaction between vortices. In fact it is rather a dy- 
namical effect of the CS term - in the limit of very small k 
vortices become very thin - the mass defect goes to zero. On the 
other hand for k = 0, in the nonlinear Schrodinger equation, the 
core width is infinite. The limit n — > is singular - one can not 
analytically continue Eq. ([!{]) to vanishing ft. The crucial effect 
of the Chern-Simons gauge field is to make the mass defect finite 
- all the fields around a vortex are well localised. 

Vortices appear to be anyons in the gauged nonlinear 
Schrodinger equation with Chern-Simons or Chern-Simons- 
Maxwell term. One might wonder if it is also possible to have 
anyons with just the Maxwell term. The answer is negative. The 
phase factor is proportional the mass missing in the vortex core 
(p is less than 1 there). However in the Maxwell- Schrodinger case 
p is also charge density and because of screening the net charge 
(mass defect) of a vortex must be zero. Addition of the Chern- 
Simons term modifies the Gauss' law and nonzero net charge 
becomes possible. 

Vortices in type II superconductors. The fact that 
Maxwcll-Schrodingcr dynamics does not admit anyons does not 
rule out anyons in real superconductors. The evolution of the 
condensate wavefunction is described by such a model but there 
are normal charges and currents in addition to those of the con- 
densate. Because of screening the total net vortex charge has to 
vanish. Such a weakened restriction does not rule out the pos- 
sibility of nonzero net superconducting charge and that is what 
really matters. Below we will derive statistical interaction be- 
tween vortices with an assumption of local charge neutrality but 
this assumption can be released in the dilute vortex gas limit. 

The time-dependent Ginzburg-Landau model is defined by the 
following Lagrangian density 



L = -iF^F^ + liWDot-cc.) 
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The Lagrangian is already in rescaled dimensionless units. We 
will consider this model as 2 + 1 dimensional and later on dis- 
cuss modifications in for real thin superconducting films. A is 
a dimensionless constant which in the considered case of type 
II superconductors is larger then 1. The covariant derivatives 
are D u ip — d a ip + iA^ip. po is an uniform positive background 
charge density which in our rescaled model is equal to +1. J n is 
the mentioned normal current minimally coupled to the gauge 
field. 



Variation with respect to Aq leads to Gauss' law 

- V 2 A° + d t (divA) =p b - + J° ■ (21) 

Now we assume local charge neutrality. In other words we assume 
the distribution of the normal current is governed by some fast 
degrees of freedom and it quckly adjusts itself so that the R.H.S. 
of the above equation vanishes locally. This assumption is far 
too strong for our needs and it is going to be released. With this 
assumption we can choose the Coulomb gauge divA = together 
with Ao = 0. We can forget now about the scalar potential and 
thus the model in a static case reduces to the Ginzburg-Landau 
theory 
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It is well known jlTj that such a model admits well localised 
vortex solutions. Field equations are 



V 2 A l 



1 



i(ip*Ditp - c.c.) 



D k D k i> = - \(i>*i> - 1)V> 



With an Ansatz 



iP(r,6) = h n (r)e m » , A» (r, 6) = - A n (r) 



(23) 



(24) 



a vortex solution with a winding number n is obtained. Both the 
scalar field and the magnetic field are exponentially localised. For 
A > 1 the penetration lenght of the magnetic field is larger then 
the coherence lenght of the scalar field. The gauge invariant com- 
bination d k ip+iA k ip is exponentially falling down with a distance 
from a vortex core. Multivortex solutions can be approximated 
by a product Ansatz the same as in Eq.(^) supplemented by an 
additive Ansatz for the gauge fields. 

Once again we can consider Berry phases picked up by the 
wave-function thanks to the Schrodinger term. This term is still 
the only term which can contribute to the Berry phase because 
the only other term with time derivatives in the Lagrangian J20| ) 
is already quadratic: ^d t A k d t A k . Thus the considerations as to 
the Berry phase in the second section can be directly applied to 
the present model with a result 
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where this time Mq is a well-defined integral 
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being a total mass defect of the superfluid directly related to the 
missing charge Q = —M . The last equivalence in Eq.(??) is a 
definition of the radius of the vortex core. In this model vortex 
is a composite of a nonzero net superfluid charge and a nonzero 
net magnetic flux. One might think that this is the origin of the 
Aharonov-Bohm effect. In fact it is just a mere coincidence. 

Another effect considered before in a microscopic setting |Q 
is a Berry phase for a single vortex ( [ll"||l2] ) 
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This Berry phase is responsible for the Magnus force. 

These are not all the terms in the effective action. Another 
term is a potential interaction derived in []l9f which is expo- 
nentially suppressed for widely separated vortices. Finally, as 
we already have mentioned, there are terms in the Lagrangian 
( |20| ) quadratic in time derivatives ^d t Akd t Ak- These terms con- 
tribute to a term in the effective Lagrangian quadratic in veloc- 
ities 



dt 2 me ffXvX v 
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One could naively think that m e f t can be found by substituting 
to T^dtAkdtAk multivortex fields with time-dependent positions 
of vortices. In fact one would have to take into account deforma- 
tions of the multivortex fields for a given trajectory up to terms 
linear in velocities. The calculations are involved but they have 
been performed for relativistic Chern-Simons vortices |l(J. We 
are going to persue this topic in a separate publication. There 
is an intriguing possibility that rn e ff can be determined by the 
parameters of the model (pf)|). 

Thus we have shown that unlike vortices in superfluid he- 
lium films those in superconducting films are anyons. In fact we 
have considered a strictly 2 + 1 dimensional model. In reality 
superconducting films are immersed in 3-space. Vortices in such 
quasi-planar systems are not so well localised. The combination 
dktp + iA^tp outside of the core does not fall down exponctially 
but rather like 4^ J2(|. In superfluid helium dkip falls like i. The 
decay rate in a superconding film is strong enough for the inte- 
gral in Eq.(p6|) to remain convergent. Thus our conclusions are 
qualitatively unchanged for thin superconducting films. 

Finally let me stress that the local charge neutrality assump- 
tion is not essential. The essential point is that there is a nonzero 
net superconducting charge associated with a vortex. This net 
charge is enclosed by the trajectory of another vortex. It does 
not matter if there are any local charge fluctuations provided 
they are localised close to the core of the vortex. 

Extension to 3+1 dimensions. So far we have considered 
only planar systems. Our goal was to make our discussion some- 
what parallel to the literature around the quantum Hall effect. 
However because our anyons arise without any help from any 
Chern-Simons term one is free to think about generalisations to 
3+1 dimensions. Of course there is a trivial embedding of pla- 
nar solutions in 3-space. There should be statistical interaction 
between parallel vortex lines. From a theoretical point of view 
statistical interaction between vortex loops would be more in- 
teresting. One could consider closed worldsheets of the loops. If 
it were not relevant one could similarly as on a plane define a 
reference loop and close any open worldsheet by continuing its 
boundaries to the chosen reference loop. 

Similarly as in the planar case one can work out the Berry 
phase for vortex loops parametrised by time t and lenght-like 
parameters 07 
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S eff = ^rJ2 dtdai x X'l) 



where I runs over the loops. " ' " means differentiation with 
respect to the appropriate parameter a. The paramctrisation by 
ct's is so chosen that the circulation of the phase is clockwise 



when one looks along the tangent to the loop. This is the form 
of the Berry phase in the limit of vanishing core thickness. The 
modulus of the phase for a closed loop worldsheet is 2tt times 
the volume enclosed by the worldsheet. Similarly as on the plane 
one has to be careful with its sign. 

Mutual statistical interactions arise when one takes into ac- 
count the finite core thickness. The general formula for the Berry 
phase in a compressible case (finite core thickness) reads [[| 
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If the core thickness were neglected one could put p = 1 every- 
where to obtain Eq.(p9|). Thus the modification to Eq.(p9|) comes 
only from vortex cores where p deviates from the vacuum value. 
The Berry phase in 3 dimensions is not a topological index as 
can be seen in the following example. 

Let us consider a pair of loops. One of them is a circle of radius 
R in the x-y plane centered at the origin 



X\(t,ai) = (R cos g\ , R sin o\ , 0) 



(31) 



The second loop is also a circle in a plane parallel to the x-y 
plane but its radius and center are time-dependent 

X2(t, 17 2) = [{R + a cos cot) cos (T2, (R + acosut) sin 02, asmujt] . 

(32) 

The worldsheet of the second loop is a torus with the first 
loop inside of it. The only contribution to the Berry phase 
comes from the motion of the second loop and is equal to 
— 2tt x ("volume of the torus" — 2it 2 RR 2 ). The second term 
is the volume of the thin vortex core of the first loop. This term 
is an analogue of the statistical interaction of vortices on a plane 
- it is the mass missing inside the closed worldsheet. This term 
is not a topological invariant because its value depends contin- 
uously on the radius R while the topology is R- independent. 
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